Cooperative Charging in a Nanocrystal Assembly Gated By Ionic Liquid 
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In order to make a densely packed assembly of undoped semiconductor nanocrystals conductive, 
it is usually gated by a room temperature ionic liquid. The ionic liquid enters the pores of the 
super-crystal assembly under the influence of an applied voltage. We study the capacitance of such 
a device as a function of the gate voltage. We show that, counter-intuitively, the capacitance of the 
system is the sum of delta-functions located at a sequence of critical gate voltages. At each critical 
voltage every nanocrystal acquires one additional electron. 



I. INTRODUCTION 



Densely packed assemblies of monodisperse semicon- 
ductor nanocrystals (NCs) are intensively studied as new 
materials with tunable optical and electronic properties. 
They potentially can be used for solar cells, light emit- 
ting diodes, transistors and photodetectors [ll, 0]- When 
NCs are not doped and, therefore, are insulating one 
can gate them electrochemically for example, by a 
room temperature ionic liquid (RTIL). Electrochemical 
gating allows one to bring several electrons per NC [1]. 
It was shown that charge transport through the assembly 
is characterized by variable range hopping, whose param- 
eters change non-monotonically with the gate voltage V 
This calls for an understanding of the charging pro- 
cess, namely a calculation of the capacitance of a densely 
packed NC assembly as a function of V, which is absent 
in the literature. 

In this paper, we are concerned with the charging 
process of a multi-layer, three-dimensional super-crystal 
(SC) of semiconductor NCs resting on a grounded metal- 
lic contact (see Fig. [1]). The SC is covered from above 
and gated by a bulk RTIL in which a second metallic 
contact is immersed. The RTIL is a good conductor, so 
that when the voltage V is applied between the top and 
bottom metallic contacts, anions from the RTIL are at- 
tracted to the top electrode and cations are attracted to 
NCs that are negatively charged by electrons. As a result, 
there are effectively two series capacitors, one formed by 
the top metal-RTIL contact and the other formed at the 
bottom RTIL-semiconductor SC contact. Both contacts 
are blocking with respect to transport of ions from the 
RTIL and electrons into the RTIL over a large window 
of \V\, so that no insulator is necessary. 

In a typical experiment the contact area of the top 
metallic electrode is made so large (not shown in Fig. [1]) 
that the latter capacitance plays no role in the total ca- 
pacitance between the two gating electrodes, since it is 
placed in series with the much smaller SC-RTIL capac- 
itance. In other words, all the external gate voltage V 
drops at the SC-RTIL capacitor, which is the object of 
our interest below. We assume that all NCs are identical 
spherical particles of internal radius a and are covered 
by a thin layer (^ 0.5 nm) of an insulator to prevent 
them from being sintered. We further assume that the 
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FIG. 1. A section of a three-dimensional face-centered cu- 
bic crystalline assembly (super-crystal) of semiconductor NCs 
(blue spheres, light gray) resting on a grounded metallic con- 
tact (yellow, very light gray). The SC is covered from above 
by a RTIL (green, gray) in which a second metallic contact 
(yellow, very light gray) is immersed. When a voltage V is 
applied between the top and bottom metallic plates, cations 
from the RTIL (small red spheres, dark gray) fill the octahe- 
dral pores of the SC to compensate the negative charges of 
electrons arriving to NCs from the bottom electrode. 



external diameter (including the insulating layer) of NCs, 
2& ~ 5 nm (a and b are shown in Fig. [2]), is large enough 
so that both cations and anions of the RTIL, which have 
diameter ~ 1 nm, can freely percolate from the bulk of 
the RTIL into the pores between spheres of the SC and 
thereby establish unique, uniform electrochemical poten- 
tials. Thus, the RTIL can be considered as an equipo- 
tential ionic conductor. It effectively serves as one of the 
electrodes of the bottom capacitor we would like to study. 
The other electrode is the SC itself. 

We assume that the electrons of all NCs in the SC have 
the same electrochemical potential as the bottom metal- 
lic electrode in spite of the fact that the semiconductor 
NCs are covered by an insulating layer. This implies 
that tunneling across the barrier between two neighbor- 
ing NCs is fast enough that electrons always have enough 
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FIG. 2. Cross section of the face-centered cubic lattice of 
NCs along the (100) plane. One can see the rock salt lattice 
of singly charged NCs (large blue spheres, gray) and cations 
(small red spheres, dark gray). The internal diameter 2a and 
external diameter 26 are shown. The arrows show a possi- 
ble choice of dipoles made from nearest neighbor cations and 
electrons. 

time to hop between NCs many times during the period 
of the AC voltage signal used for capacitance measure- 
ment. In this sense, the bottom metallic electrode and 
the SC together constitute the equipotential electrode of 
the bottom capacitor. Again, the RTIL is the second 
electrode. On the other hand, we also assume that the 
electron tunneling barrier between neighboring NCs is 
large enough to allow only an integer number of electrons 
in each NC. This requires the tunneling matrix element 
to be smaller than any residual disorder distinguishing 
between neighboring NCs. 



II. NONTRIVIAL CHARACTER OF 
CHARGING OF SC-RTIL CAPACITOR 

Let us now consider charging of the SC-RTIL capaci- 
tor. We assume that semiconductor NCs are not doped, 
so that some gate potential is necessary to align the 
Fermi level of electrons with the conduction band of the 
semiconductor. In this paper we calculate the gate volt- 
age V from this reference point. Our goal is to calcu- 
late the three-dimensional concentration of electrons n 
as a function of V as well as the volumetric differen- 
tial capacitance (capacitance per unit volume) of the SC: 
C{y) — edn/dV. Of course, the validity of this concept 
requires that, e.g., the capacitance of the SC be propor- 
tional to the number of NC layers in the NC film. 

It is natural to think that the charging of all NCs of 
the SC assembly happens sequentially with growing V, 
i.e., NC by NC, because as in any capacitor an electron 
just added to a particular NC repels other electrons and 
prevents them from entering neighboring dots. This in- 
deed would be true if we were dealing with a single layer 
of NC charged by parallel, planar metallic gate hanging 
above it and separated by a thin insulator (Fig. [3]). In 



FIG. 3. An example of a more conventional capacitor. One 
layer of NCs (blue spheres, gray) is charged by parallel plane 
metallic gate (yellow, light gray). Electrons within the NCs 
make dipoles with counter charges (arrows) in the top gate. 
All dipoles are collinear and reside in a single plane perpendic- 
ular to the two parallel metal plates. The repulsion between 
dipoles leads to monotonous growth of n{V). 



such a system each electron arriving to a NC makes a 
dipole with its image in the top gate. All dipoles reside 
in a single plane perpendicular to the two metal plates. 
They, therefore, repel each other and at small concentra- 
tions form a correlated liquid of dipoles or dipole Wigner 
like crystal, which provides resistance to charging of the 
NC layer. The capacitance of such a system was calcu- 
lated in Ref. [5] and at finite temperature it is everywhere 
positive and finite, so that n(V) is a monotonically grow- 
ing function. In this sense, the system shown in Fig. [3] 
is qualitatively similar to a conventional plane capaci- 
tor made from two metallic parallel electrodes, the en- 
ergy of which may be calculated as the result of pairwise 
repulsive interactions between elementary, continuously 
distributed dipoles created by positive and negative in- 
finitesimal charges on two opposite capacitor plates. 

In this paper, however, we predict that in a SC-RTIL 
capacitor the electron density n{V), or, equivalently, the 
filling factor of NCs i'{V) = n{V)/N, shows alternat- 
ing long plateaus and vertical steps. At each step n{V) 
changes by the concentration of NCs N or, in other 
words, each NC of the SC acquires at the same V ex- 
actly one additional electron. This means that C{V) is 
the sum of a small number of delta-functions located at 
special voltages Vk , at which the fc-th electron enters each 
dot. Schematic plots of the filling factor i/ and capaci- 
tance C of the SC as a function of gate voltage V are 
shown in Fig. 01 Together, the whole SC behaves like a 
single quantum dot in spite of the interaction between 
NCs. Indeed, in a single NC or a quantum dot the fc-th 
electron is admitted at the critical voltage Vk = ke/C, 
where C is the gate-dot capacitance. The voltage Vk sat- 
isfies eVk — Ey=k ~ -Et)=fe-i, where E^=k is the energy 
per NC in a SC with filling factor v = k. 

Throughout this paper, we use the term "cooperativ- 
ity" to refer to the fact that all NCs of the SC simulta- 
neously acquire an extra electron at a particular voltage 
as a consequence of an attractive net interaction between 
electron-cation dipoles. Thus, our goal in this paper is 
to prove and understand cooperativity. We concentrate 
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dral pores see Ref. [e"].) There are two tetrahedral pores 
and one octahedral pore per NC. The octahedral pore 
can admit a sphere with radius r = (\/2 — 1)6 ~ 0.415, 
while a tetrahedral one can admit a sphere with radius 
r = (a/6/2 - 1)6 ~ 0.225. Let us first assume that NCs 
are relatively small and cations can only fit into the oc- 
tahedral pores. In this case, cations can fill up all the 
octahedral pores of the SC in order to make the assem- 
bly overall neutral (Fig. [1]). Each cation is assumed to 
be localized in the center of an octahedral pore. A single 
electron in a NC is confined by the NC wall and is also 
affected by the Coulomb potentials of cations and neg- 
ative charges of other NCs. We assume that the radius 
a of the semiconductor NC is much smaller than the ef- 
fective Bohr radius as = h^e/me^ of the semiconductor. 
Here m is the effective electron mass in the semiconduc- 
tor. In this case, in the first approximation the kinetic 
energy dominates the Coulomb energy and we may use 
the unperturbed electron ground state wave function 



FIG. 4. Schematic plots of the filling factor of NCs v{V) = 
n{V)/N and the SC-RTIL capacitance as a function of volt- 
age V . i^{V) shows alternating long plateaus where it is an 
integer and vertical steps, where i^iV) changes by unity. The 
capacitance C is the sum of delta-functions located at special 
voltages Vk, at which k-ih electron enters each dot. 



mostly on the case of the transition between i> — and 
= 1. 

In Sec. Illli we study a simple model of semiconductor 
NCs with large Bohr radius as, i.e., weak Coulomb inter- 
action, and demonstrate the existence of cooperativity. 
In Sec. IIVI metallic NCs are discussed and are proven 
to show cooperativity of electrons as well. Appendices 
A, B, and C deal with computational and mathematical 
details. 



III. SEMICONDUCTOR NC ASSEMBLY 

In this section, we focus on a densely packed semicon- 
ductor NC assembly. Although in a real system the NCs 
usually have a different dielectric constant e than that of 
the RTIL, as a pedagogical example to illustrate cooper- 
ativity, e is for the moment taken to be uniform every- 
where. More realistic cases with non-uniform dielectric 
constant will be discussed later in this paper. We show 
below that the counter-intuitive behavior of the SC-RTIL 
capacitor arises from the combination of two conditions: 
1) the electrodes of the SC-RTIL capacitor interpenetrate 
each other, and 2) charges of both signs are discrete and 
localized. 

Let us assume that we deal with a 3D face-centered 
cubic SC of densely packed NCs (Fig. [T|). Within an 
fee lattice, there are two types of pores, octahedral and 
tetrahedral. (For visualization of octahedral and tetrahe- 
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and the corresponding ground state energy per NC 
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For a calculation of the main term of the electrostatic 
energy of the system, one can think that each electron is 
located in the center of the NC. Octahedral pores form 
an fee lattice that is shifted in space relative to the fee 
lattice on which the NCs reside. These two alternating 
fee lattices, the one of singly charged NCs and the one of 
cations, together form a rock-salt lattice of point charges. 
Because of complications of showing a three-dimensional 
lattice, we illustrate a mutually penetrating SC of NCs 
and the crystal of cations for the case of a 2D square 
rock-salt lattice SC at 1/ = 1 (Fig. [5]). 

One may divide the rock salt crystal into neutral pairs 
of a cation and an electron in neighboring NC and think 
that these pairs enter sequentially as the voltage V in- 
creases. However, the following consideration shows that 
such pairs enter simultaneously at some threshold voltage 
rather than sequentially. Let us first calculate the volt- 
age Vg at which one such neutral pair enters an initially 
empty lattice. Bringing one electron into a NC costs 
quantum mechanical energy Eq, plus some energy Eion 
to bring a cation into a pore of the empty system from 
the bulk of the RTIL. Generally, Eion is determined by 
the chemical potential of ions in the bulk RTIL, which is 
of the order of the nearest neighbor interaction between 
ions in the bulk RTIL. In addition, one should include the 
Coulomb interaction between the electron and the cation: 
-e^/e^/2h. Thus, the total energy associated with bring- 
ing one electron-cation pair into the lattice is 
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At first glance it seems natural to think that when the 
voltage V is increased to eVs = Eg the first electron- 
cation pair enters an empty SC. However, this is incor- 
rect, as can be seen by considering the total energy of the 
V = \ state, at which each NC of the SC is occupied by 
a single electron that is compensated by a neighboring 
cation. In this way the whole rock-salt lattice of singly 
charged NCs and cations is formed. The electrostatic 
energy per pair in this lattice is actually 



Em, =-Mi- 



(4) 



where Mi ~ 1.75 is the Madelung constant for a rock-salt 
lattice f7|. Accordingly, the total energy per electron- 
cation pair in the filled SC lattice is 



Eu=l = Eq+EMi+Eion ~ ^—1.75 — —+Eu 
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Comparing Eq. ([5]) to Eq. Q we see that for the fully 
filled lattice the energy per pair is smaller by 0.51e^/e&. 
Therefore, when the voltage V is tuned up, no electron- 
cation pairs enter the empty SC until the critical voltage 
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at which N such pairs enter collectively. Thus, charging 
of the SC happens as a first order phase transition. As 
a consequence, the capacitance has a 5-function peak at 
V ^Vi (Fig.m). 

One can also use the language of dipole moments 
formed by nearest neighboring electron-cation pairs in 
the SC-RTIL capacitor to illustrate the mechanism of 
this first order transition. In the rock salt arrangement of 
electrons and cations, apparently the repulsion between 
neighboring, aligned dipoles is overcome by the attraction 
between anti-aligned dipoles (Fig. [5]). This happens be- 
cause the vector connecting the centers of nearest neigh- 
boring dipoles forms an angle of 7r/4 with respect to the 
dipole moment. In this sense, the SC-RTIL capacitor in 
the form of a rock-salt lattice of electrons and cations is 
dramatically different from a more conventional capaci- 
tor where all dipoles repel each other (Fig. [3]). It is the 
dipole-dipole attraction that leads to cooperativity in the 
charging of the SC-RTIL capacitor. 

For the SC of small NCs discussed above, the cations 
can fit into the octahedral pores only. However, for a SC 
of larger NCs, where the pore size is thereby larger, the 
lattice can accommodate cations in its tetrahedral pores 
as well and a second delta peak of capacitance might oc- 
cur. We show below that, for the case of cations residing 
in the tetrahedral pores, the system can exhibit cooper- 
ativity aX V = 1 and v = 2 only. 

Recall that there are two tetrahedral pores per NC 
within an fee lattice. For singly charged NCs, the cations 
fill half (one sublattice) of the tetrahedral pores instead 
of the octahedral ones. The corresponding lattice of 



charges resembles the zinc blende lattice such as ZnS (6[ , 
where each NC takes the place of S^~ and cations occupy 
Zn^+ sites. To bring a single electron-cation pair into the 
empty SC, it costs energy 
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where the first term is the electron ground state energy 
Eq, the second term represents the Coulomb interaction 
between the electron and its nearest neighboring cation 
(•\/6/26 is the distance between them), and Eion is the 
chemical potential of cations in the bulk RTIL. For the 
fully filled SC where each NC is occupied by an electron, 
the total energy per electron-cation pair is 
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where L64 is the Madelung constant for a zinc blende 
lattice [7|. 

Comparing Eq. ([7]) with Eq. ([5]) we find that energy per 
pair is smaller by 0.52e^/e6, which means that the first 
electron enters each NC simultaneously at eVi = E^^i, 
and that at Vi, the capacitance C{Vi) shows a delta peak 
(Fig. 13]). Comparing Eq. ^ with Eq. we see that a 
cation indeed prefers a tetrahedral pore to an octahedral 
one. The reason is that in the former pore the cation is 
closer to its neighboring electron than in the latter. 

When a second electron enters each NC, it occupies the 
same Is state but with opposite spin. To make the assem- 
bly overall neutral, cations from the bulk RTIL fill the 
other half of the tetrahedral pores so that there are two 
cations per NC. This lattice of point charges is equivalent 
to the fluorite CaF2 lattice. Each NC with two electrons 
takes the place of Ca^"*" and cations occupy F~ sites. 

To see whether the system exhibits cooperativity at 
v = 2, one should compare the energy of a single addi- 
tional electron-cation pair with the energy per pair of the 
state where N additional electrons enter the zinc blende 
lattice collectively. When only one additional electron en- 
ters, a cation resides in the originally empty neighboring 
tetrahedral pore to compensate for it. The energy cost 
to bring this pair into the zinc blende lattice is therefore 

c2 
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where Er is the repulsion between two electrons within 
a NC, Ec + Er is the electrostatic interaction energy be- 
tween the zinc blende lattice of singly charged NCs and 
cations in first sublattice of tetrahedral pores and the 
newly-inserted nearest neighboring electron-cation pair, 
the last term is the Coulomb interaction between the 
pair, and Eg and -Eion are defined in the same way as 
above. In contrast, when N additional pairs enter the 
SC, they form a second zinc blende lattice and the total 
energy per such pair is as follows: 

„2 
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By comparing Eq. ([9]) with Eq. ^TU\\ we see that again the 
energy per additional pair is lower than that of a single 
additional pair. This means that cooperativity exists at 
1^ = 2. 

Since Ec is unknown, to obtain the voltage V2 at which 
a second electron enters each NC simultaneoulsy, let us 
first calculate the total energy per pair of the SC a.t 1/ = 2. 
This energy is the sum of the quantum mechanical en- 
ergy 2Eq, the direct Coulomb repulsion Er between two 
electrons within a NC, the electrostatic energy per pair 
of two electrons and two cations Ems^ and the chemical 
potential of the cations 2Eion'- 

Eu=2 — '^Eq + Er + Em3 + 2Eion (H) 

= T + 1.79 5.04—^ + 2E,on. 

ma'^ ea £ 2 ^ 

Here, 5.04 is the Madelung constant for a fluorite lattice 
0. Therefore, 

eV2 = E,=2 - E,=i - + 1.79— - 3.4—^ + E,on, 

Ima'^ ea e—b 

(12) 

where a second delta-function of capacitance is located 
at potential V2. 

For higher filling factor like ly — 3, since both Is states 
are already occupied, a third electron entering NC will 
be in the Ip state. This quantum mechanical energy is 
greater than the Coulomb interaction of electron-cation 
pairs and thereby causes a large gap between V2 and V3. 
This gap also increases because one can expect larger 
energy necessary to bring a third cation per NC into the 
assembly. To take into account this effect one should 
change the model and go beyond one cation per pore or 
consider SC lattice swelling. We have not tried to explore 
such more complicated models yet. 

IV. METALLIC NANOCRYSTAL ASSEMBLY 

Throughout Sec. Ill, the dielectric constant e was re- 
garded as uniform over the whole lattice. This is an over- 
simplification of a real system, since semiconductor NCs 
usually have e larger than that of a RTIL. As a conse- 
quence, one can expect that cations create image charges 
in surrounding NCs, which complicates the calculation of 
energy as a function of cation concentration. To under- 
stand whether image charges can qualitatively affect co- 
operativity, in this section we study metallic NCs, which 
corresponds to the case where NCs have infinite dielec- 
tric constant, where the effect of images should be even 
stronger. Specifically, we consider a SC of electrically iso- 
lated spherical metallic NCs with external radius b (see 
Fig.[2|). 

Suppose that the insulating film covering NCs makes 
the overlap between electron wavefunctions of neighbor- 
ing metallic NCs so weak that the charge on each NC is 
quantized and the SC conducts only by electron hopping. 



Metallic NCs are assumed to form a face-centered cubic 
lattice and cations reside in the octahedral pores (Fig.[T|). 
We show below that for this system the cooperativity of 
electrons still exists for the case of metal NCs in the vicin- 
ity of filling factor ly = 1. This conclusion suggests that a 
SC of semiconductor NCs (with large but finite dielectric 
constant) gated by RTIL should also show cooperativity 
near filling factor v = I. 

Let us first calculate the electrostatic energy of the sys- 
tem occupied by a single pair of a cation and an electron 
in its adjacent NC. By classical electrostatics we know 
that the cation will induce images in all metal NCs of 
the SC [1]. These image charges produce new images in 
every other metal sphere, and so forth. Thus we have 
an infinite series of image charges in each sphere. The 
problem of image charge summation has previously been 
discussed for the case of a point charge in the vicinity of 
two conducting spheres [9^; here we discuss the problem 
of an infinite lattice of conducting spheres. 

The electrostatic energy cost to bring a single electron- 
cation pair into the empty lattice is determined by the 
total interaction energy of the cation with all the in- 
duced images. It should be noted that for a cation in 
the center of an octahedral pore the charges created in 
the cation's six nearest neighboring NCs are dominat- 
ing. We computed the electrostatic interaction energy 
between one cation and the image charges induced in its 
six nearest neighboring NCs. To solve this problem nu- 
merically, we construct an infinite series of image charges 
via a recursion process (see Appendix A). The compu- 
tational result shows that the electrostatic energy of an 
electron-cation pair in the empty lattice is approximately 
Eg = (— 1.06± 0.1)e^/b, where b is the external radius of 
NC. 

We now turn to the calculation of the energy per pair 
in a fee lattice of singly charged metal NCs with cations 
in all the octahedral pores. Since each cation creates im- 
ages in every metal sphere of the SC, the image charge 
summation becomes very difficult to carry out. In order 
to solve this problem we adopt the following numerical 
strategy. We first replace each cation by a metallic sphere 
with vanishingly small radius Oc and charge +e. We then 
numerically solve the Laplace equation at a given V by 
the Gauss relaxation method [Ifll- From this numeric 
solution one can determine the charging energy, namely 
the total energy of the system per electron-cation pair 
(see Appendix B for more details). To obtain the inter- 
action energy per such pair we need to subtract the self- 
energy of a metal sphere e^/2b and of a cation e^/2ac 
in isolation from the charging. Numerical results give 
E,,^i = (-1.48 ± 0.01)e76 at < b. 

Comparing Es with E^^i, we see that the electrostatic 
energy per pair in the filled lattice is lower than that of a 
single pair, which means that there is indeed cooperativ- 
ity of electrons in a filled metallic SC. This result can be 
directly applied to a SC of heavily doped semiconductor 
NCs, which behave as metallic dots. On the other hand, 
it also implies that for a SC of undoped semiconductors 
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the large dielectric constant of the semiconductor does 
not suppress cooperativity in the charging process of a 
SC gated by RTIL. 

In order to test the generality of these conclusions, we 
examined the case in which metallic NCs form a simple 
cubic lattice instead of an fee lattice. For a simple cu- 
bic lattice, each unit cell contains one cation residing in 
the center and eight nearest neighboring NCs occupying 
the corners. The methods described in Appendix A and 
B are employed to numerically compute the electrostatic 
energy of a single electron-cation pair and the energy 
per pair in a filled lattice. As expected, the simple cu- 
bic lattice of metallic NCs exhibits cooperativity as well: 
the electrostatic energy per pair of a fully-filled lattice 
is approximately — 1.22e^/&, lower than that of the sys- 
tem occupied by a single pair, which is about — 0.63e^/6. 
This again corroborates our prediction that cooperativ- 
ity is a general property of electrons in nanocrystal SCs 
rather than a phenomenon applicable only to a few spe- 
cial cases. A more general proof of cooperativity of SC 
of metallic NCs is given in Appendix C. 

Until now we have not discussed the difference between 
the charge distributions in metallic NCs and semicon- 
ductor NCs with large dielectric constant e. For metallic 
NCs, all charges are distributed on the surface, while 
for semiconductor NCs, the injected electron is confined 
within the NC with a ground state wave function given 
by Eq. ([ij. To see whether this difference affects coop- 
erativity, let us first look at a singly charged NC and its 
adjacent cation. For the case of a metallic NC, the sur- 
face charge of the NC can be seen as a superposition of 
two parts: (1) charge that is induced by the cation in a 
grounded NC, and (2) an electron charge — e that is uni- 
formly distributed on the surface. For a semiconductor 
NC, the total charge is also composed of two parts, but 
the second part is an electron with the spherically sym- 
metric ground state wave function given by Eq. ([T]). In 
the semiconductor case one should pay an additional self 
energy Eas to assemble the electron distribution implied 
by Eq. ([IJ from the uniform surface charge which has to- 
tal charge — e. Eas is the same for each NC and therefore 
is independent of whether the lattice is empty or filled. 
When considering a fully filled lattice of semiconductor 
NCs, since every NC has such an energy difference due to 
charge redistribution, the critical voltage at which coop- 
erativity shows up is in fact shifted by Eas/e compared 
to the result of metallic NCs. Apparently this shift of 
threshold voltage does not have any influence on coop- 
erativity. Therefore, our analogy between metallic NCs 
and semiconductor NCs with large dielectric constant is 
conclusive. 



V. CONCLUSION 

In summary, we have shown that the capacitance of 
an ideal supercrystal of semiconductor NCs gated by a 
RTIL is the sum of delta-functions located at critical gate 



voltages. At each such critical voltage every NC aquires 
an additional electron. This was shown assuming that 
cations arrive in the middle of the pores of the super- 
crystal and arrange themselves in a crystalline lattice. 
To our mind this is an interesting theoretical model of 
a macroscopic system of interacting nanoparticles which 
demonstrates the capacitance behavior expected for an 
ensemble of non-interacting quantum dots. 

In our theory we ignored the effects of finite tempera- 
ture and disorder on the crystalline lattice of cations and 
charged NCs. The thermal energy kgTji corresponding 
to room temperature Tfj should be compared with the 
gain which a nearest neighbor pair of a charged NC and 
a cation gets when cations and electrons enter cooper- 
atively. This gain W ~ 0.5e^/eb, as we saw above for 
the system with uniform dielectric constant. For a semi- 
conductor NC with dielectric constant much larger than 
the dielectric constant '-^ 2 of the RTIL, the energy 
W ~ 0.5e^/eib ^ bkBTn for b = 2.5 nm. This means 
that temperature plays a secondary role in the thermo- 
dynamics of our system and the capacitance still should 
have delta function peaks. 

Now let us discuss the role of a spatial disorder. One of 
the sources of such disorder can be size fluctuations 5a of 
the NC radius a. In experiment the relative fluctuations 
5a/ a are typically around 7% These fluctuations of 
course affect the structure of the SC and all Coulomb en- 
ergies. But for relatively small NCs the most important 
role of size fluctuations is the variation of the large quan- 
tum mechanical energy of electrons within the NC. For 
a = 2 nm and relative size fluctuation — 7%, for example, 
we get 

5Eq = Eq'^ = O.UEq ~ Q.leV ^-ikBTR, (13) 

Such fluctuations are comparable with the energy gain 
associated with cooperativity and may destroy cooper- 
ativity effects and smear out the delta peaks of capac- 
itance predicted above. Instead, if smaller fluctuations 
of size of order of 4% [ll| are achievable one can count 
on the observation of delta function peaks in the capac- 
itance. However, we are not aware of any experimental 
demonstration of such a capacitance. 

The most advanced study of charging of a SC of al- 
most identical semiconductor NCs Q focuses on the de- 
pendence of conductivity of the SC, (t(T), on the gate 
voltage V. The authors discovered that the conduc- 
tivity follows the variable range hopping (VRH) law 
cr(T) cx exp[— (To/T)*], where the power s oscillates with 
growing V from s = 1/4 (Mott VRH [11]) to s = 1/2 
(Efros-Shklovskii VRH 1131 ). (Such periodic oscillations 
were predicted in Ref. |14| for periodic granular system 
of metallic dots doped by donors that are randomly po- 
sitioned in the spaces between the dots.) Apparently, 
observation of variable range conductivity in the samples 
of Ref. [4] points to the important role of disorder. This 
disorder may be a random distribution in sizes, which 
leads to fluctuations in the quantum energy, as discussed 
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above. One can also imagine that in a system with nar- 
row pores cations never become completely equilibrated 
to their crystalline ground state arrangement. These is- 
sues should all be clarified by future experiments, where 
capacitance data should be analyzed together with trans- 
port one. 
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Appendix A: Image Charge Summation 

In this Appendix we show how an infinite series of im- 
age charges is formed in each metal sphere of the SC sur- 
rounding a point charge. Consider a cation in the octahe- 
dral pore surrounded by six metal spheres (metal NCs). 
Let us take the singly charged NC adjacent to the cation 
as an example. The cation first creates an image charge 
—e' = —eh/^/2b = — e/\/2 {\/2h is the nearest neighbor 
distance) lying on the ray from the center of a neighbor- 
ing NC to the cation, as described in Ref. [1]. To main- 
tain the total charge — e in the metal sphere and equipo- 
tential on the surface, a second image e" = —e + e/\/2 
has to be placed at the center Q. These two charges 
induce images respectively inside all the other spheres. 
The newly induced images in each neutral NC will again 
create new images in the singly charged NC, and so forth. 
Meanwhile, the same process is repeated for all the other 
NCs, except that the total charge in each of them is zero 
instead of — e. In this way an infinite series of image 
charges is formed inside every metal sphere. 

The above image generating process can be written as 
a simple recursive list of rules: 

1) Create the image charges — e', — (e — e') in the singly 
charged metal sphere. 

2) Create the image charges — e', +e' in the other five 
neutral spheres. 

3) Create images within every sphere for all image 
charges inside the other spheres. 

4) Add an additional charge to the center of ev- 
ery sphere to ensure that it has net charge — e' (singly 
charged NC) or (neutral NC). 

5) Repeat steps 3) - 4). 

The total interaction energy of the cation with the sur- 
rounding NCs can be found by summing up the Coulomb 
interaction energy of the cation with all the induced im- 
ages. The total energy of the system occupied by a single 
pair is 1/2 times this value. (Note that the factor 1/2 for 
the interaction energy arises because of the image charge 

So 

It should be noted that when a large number of metal 
spheres are included in the region surrounding the point 
charge, the above recursion process does not necessarily 
converge as a function of iteration number. (Since the 
number of image charges created during each recursion 



step grows exponentially as a function of recursion num- 
ber, and since these image charges are predominantly of 
the same sign during each step, the calculated energy 
oscillates as a function of iteration number and does not 
necessarily converge.) This problem can be circumvented 
if one introduces an artificial damping of the magnitude 
of the image charges created during step 3), such that 
image charges created during later recursion steps have 
progressively smaller magnitude. The results for electro- 
static energy presented in Sec. IIVI correspond to a damp- 
ing of the image charges by a factor 77^ , where K is the 
iteration number and 77 is some positive number satisfy- 
ing 77 < 1 and 1 — 77 ^ 1 . For the six neighboring spheres 
of the fee lattice, as considered in Sec. IIVI this damping 
is sufficient to produce convergence when 77 < 0.95. The 
convergent result for energy, Es = (—1.06 ± 0.1)e^/6, is 
seen to be independent of the choice of r/ to within the 
stated accuracy over the range 0.65 < 77 < 0.95. 

Appendix B: Gaussian Relaxation Method 

In the Gaussian relaxation method, described briefly 
in Sec. IIVI we replace each cation by a metallic sphere 
with small radius ac (usually ~ 15% of a) and charge 
+e. This charge is supported by an unknown voltage 
Ve that is applied between the cations and the metal 
NCs. That is, all cations are held at Ve while all metal 
spheres are held at zero potential. Since initially the 
value of Ve that corresponds to charges ±e is unknown, an 
arbitrary voltage V (which might not correspond to ±e) 
is assigned to each cation. We then numerically solve the 
Laplace equation at V and obtain the electric potential 
distribution in the space outside small and large spheres 
using the Gauss relaxation method. From this numeric 
solution one can find cation's charge q that corresponds 
to V, and the system's capacitance C = q/V. The total 
energy of the system per pair is essentially the charging 
energy e^/2C. Note that it is not necessary to solve Ve 
since C is linear, i.e. C — q/V = e/Ve- This procedure 
can be summarized by the following list of steps: 

1) Create a 3D lattice of grid points filling the unit cell 
of the periodic SC lattice. 

2) Enforce the boundary conditions that all lattice 
points inside a metal sphere have potential $ = and 
all points inside a cation have $ = V. 

3) Solve for the potential at all other grid points using 
the Gauss relaxation algorithm. Enforce periodic bound- 
aries at the unit cell edges. 

4) Measure the charge q of the "capacitor" comprised 
of the small (cation) spheres and the larger (NC) spheres 
by constructing a Gaussian surface around one of the 
cations and calculating § E ■ dA = Airq. The capacitance 
is C = q/V . From C we calculate the total energy per 
pair eV2C. 

Appendix C: Proof of Cooperativity for A Metal- 
lic SC. 

In the Sec. IIVI we discussed a densely packed metallic 
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SC that shows cooperativity at v = 1. In that case, the 
Gauss relaxation method is employed to calculate the en- 
ergy per electron-cation pair of a fully-occupied lattice, 
which is proven to be lower than that of a single pair 
in the system under the approximation that only the six 
nearest neighboring NCs interact with the cation. One 
may wonder whether cooperativity is a general property 
regardless of this approximation. To answer this ques- 
tion, in this section we would like to propose a proof that 
shows the existence of cooperativity for face-centered and 
cubic SC lattices of metallic NCs with large radius b hav- 
ing small metallic spheres in their pores (Fig.[T|). We sup- 
pose that the whole system is neutral and all charges are 
discrete, and in the ground state of the fully filled lattice 
all small spheres have charge +e and all large spheres 
have charge — e. Our goal is to prove in a general way 
that the energy per pair of a fully-occupied lattice, Ufun 
is lower than the energy of the system when it is occupied 
by a single +/— pair, Us- The case of point charges next 
to metal spheres can be recovered by making radius of 
the small spheres vanishingly small. 

In general, the total electrostatic energy [/ of a system 
of charged, conducting spheres is 



U 



(14) 



where the indices i and j label all spheres in the system, 
Qi is the charge of the sphere i, and C^j^ is the inverse of 
the matrix of electrostatic induction (this matrix has 
the property Cij = Cji) This expression for energy 
includes all binding energies as well as the self-energies 
of all charges. 

Consider first the case where only a single +/— pair 
is charged: say that a small sphere (cation) i ~ 1 has 
charge qi — e and a neighboring large sphere (metal NC) 
i = 2 has charge q2 = — e, while all other spheres have 
Qi ~ 0. Then the double sum in Eq. ([T4|) is reduced to 
only four terms, and the total energy Ug becomes 



22 



(15) 



Now consider the case where the lattice is entirely 
filled: all the small spheres (cations) in the lattice have 
q = e and all the large spheres (NCs) have q — — e. 
Suppose, further, that the lattice is in its lowest energy 
configuration. This implies that the system must be sta- 
ble with respect to exchanging the charges of any two 
spheres. That is, if the sign of the charge at site i and 
site j are simultaneously inverted, then the system en- 
ergy must increase or be unchanged. 

Suppose, without loss of generality, that in the ground 
state of the filled lattice the site i = 1 has charge +e and 



the site i = 2 has charge — e. From Eq. (|14p . one can 
therefore write the system's initial energy as 



U^n = 2 [e'CiV + e'C^; - 2e'C^^'] (16) 



2 E - 



J5^1,2j#l,2 



If the charges in site i = 1 and i — 2 are exchanged, 
then the energy becomes 

Uf^n = I [e'CiY + e^C^-; - 2e^C^^'] (17) 



2 {eq.C^l ^ eq,C^/) 



E E 

■1^1,2 j^l,2 



The stability of the ground state requires that C//i„ > 
Uin- By comparing Um and Ufm, one can see that this 
condition implies 



E (eqjC^^ - eq,C^l) < 0. 



(18) 



When the lattice is fully filled, the energy per pair U{ui\ 
can be found by taking the sum in Eq. (jl4p over i = 1, 2 
only, since i — 1,2 labels a single nearest- neighbor pair 
that is equivalent to all other pairs: 



"full = ^ E E'^y^^** 



i=l,2 j 



(19) 



By substituting Eqs. p3)) and into Eq. p^ . we can 
conclude that 



Ufull < Us 



(20) 



so that the lattice is proven to exhibit cooperativity. For 
the case of a metallic SC with insulating cationic charges, 
one can simply take the limit where the radius of small 
(positive) spheres is vanishingly small, and one recovers 
the case of point charges next to metallic spheres. 
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